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The Casimir effect in compact hyperbolic Universes was numerically obtained 
in previous publications. In this talk, I expose these results. 



1 Introduction 

> 

^sO ' The Universe today is very close to homogeneous. Also, no preferred direction 

OO . has been discovered, so that the most important 3-Manifolds in connection with 

| cosmology are S 3 , H 3 and R 3 . Anyway, it must be mentioned the following 

. interesting results, that isotropisation can be achieved by particle production 

■^j- | in the primordial Universe pj Gj or by the existence of a classical cosmological 

constant |3j. For the moment, let us simply suppose that the Universe really is 
""q I homogeneous and isotropic, which is in good agreement with the observations. 

U" 1 It is well known that R x S 3 , R x H 3 and R x R 3 are solutions of Einstein's 

equations (EQ). It is irrelevant for the EQ, if the spatial sections have a non 
bJ[). trivial topology modelled on them, or not. From recent determinations of the 

WMAP mission there is an indication for a spherical Universe @|, but still not 
capable of ruling out the other two cases. The low multipoles of the CMBR, are 
the ones more sensible to global properties of the space. In a very interesting 
result, showed that the low multipoles of the CMBR are better fitted, if the 
spatial sections of the Universe are the Poincare dodecahedral space, instead of 
S 3 (see also ©). 

There are many effects that arise if space is compact. The first astrophysical 
limits on the topology of the universe were obtained for a 3-torus T 3 . Accor- 
dance with the homogeneity of the CMBR puts a lower limit on the size of 
the fundamental cell, about 3000 Mpc, which is a cube in the cases of and 
jH]- Later on, it was shown that this result is very sensitive to the type of the 
compactifications of the spatial sections. For a universe with spatial sections 
T 2 x R, the fundamental cell's size is about 1/10 of the horizon, and is compati- 
ble with the homogeneity of the cosmic microwave background radiation CMBR 
jH]. In compact universes, the pair separation histogram would present spikes 
for characteristic distances. At first it was thought that this technique, known 

* muller@fis.unb.br 



1 



as the crystallographic method, was able to totally determine the topology of the 
universe ^Hj- It turned out that the crystallographic method only applied when 
the holonomy group contained at least a Clifford translation, i.e. a translation 
which moves all the points by the same distance pi] and ^2]. Generalizations 
of the crystallographic method were proposed, for example in 

Also in compact universes the light front of the CMBR interacts with itself 
producing circles in its sky pattern |14j . 

Recent results have called our attention to the possibility that methods based 
on multiple images will prove not to be efficient ^S]. The reasoning is that, ac- 
cording to observations, the curvature is very small, so the fundamental regions 
are so big that there has not been time enough for the formation of ghost im- 
ages. The result is that for low curvature universes such as ours, only compact 
universes with the smallest volumes could be detected by pattern repetitions. 

The knowledge of the spectrum of the Laplacian for compact spaces have 
some motivations. One of them is the decomposition of the CMBR in right set 
of modes, if space is non trivial. It has been studied numerically for compact 
hyperbolic space by |16U17j . And more recently QB] f° r compact spherical space. 

Casimir effect also occurs in compact spaces. In the most simple cases it can 
be obtained by analytic continuation of the zeta function. This procedure is 
known as the zeta function regularization, see for example ^Sj, [2fl| - Very elegant 
analytical results can be obtained for the Casimir effect in hyperbolic space 
|21l I22j . The method is indirect, as it makes use of the Selberg trace formula 
to obtain the zeta function, irrespective of the knowledge of the spectrum. For 
a more comprehensive introduction to the Selberg trace formula see [2H] . There 
is a recent result for the Casimir energy in shperical space [2i| . 

In this work I shall expose a different technique which was previously used 
by us GUHIIEZI. 

The paper is organized as follows. In sectional I review the simplest mani- 
festation of the Casimir effect of topological nature, and reproduce some basic 
properties of zeta functions. The static solution of EQ, the classical energy 
momentum tensor and the method of point splitting is outlined in section EJ 
In section 31 the Feynman propagator and the Hadamard function is obtained. 
And in section the Casimir energy density is numerically obtained in some 
Universes. Natural units are used, h = G = c = 1, except in section 



2 Topological Casimir Effect 

It was discovered by Casimir in 1948 [2H| that two uncharged parallel metallic 
plates in vacuum should be attracted by the force 

1 240 a 4 ' 

where S is the area of the plates and a is the distance between them. Soon 
after, it was verified experimentally [2S| - It is universal since it does not contain 
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the coupling constant of the field. Also, it is a non local effect, in the sense that 
it depends on the boundary conditions. 

Casimir effect also occurs when the topology of space is non trivial. In the 
following, let us reproduce the result for a massless, uncharged scalar field in 
the two-dimensional space Rx S 1 , for more details see [HOI GDI EH- In this case 
the spectrum of the Klein-Gordon operator reads 

lu — k = 0, k n — — 7i, n = 0, ±1, ±2, ... 
a 

where a is the volume of the manifold. The total energy is 

(b\H\b)=J2"k (n k + l 

k ^ 

where \b) is an arbitrary state of the field, with particles in mode k. The 
vacuum expectation value follows from considering — 0, and the Casimir 
energy reads 



(°i#i°> = 



2^ k 

k 



E c = ^£ 

n — ' 



a 

n=l 



As it stands, the last summation of course diverges. Nevertheless, a meaning can 
be given to it through the mathematically well known Riemann zeta function. 
Remind that the Riemann zeta function is defined for Re(s) > 1 



oo 



n—l 

This function can be analytically continued giving a finite result PI] 

The Casimir energy for Rx S 1 , then follows 

6a 

The technique used above is the famous zeta function renormalization. Appar- 
ently, Zeta functions where used in the physics literature for the first time in 
|.'S2ll.'i.'!| . In the following I shall expose the well known relation between the zeta 
function and the trace of the heat kernel. Suppose a Laplace type operator A 
acting on a compact 3-space, with a spectrum Ai, A2, ... 
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where <j> n is some complete orthonormal set of modes. The Euler infinite integral 
for the gamma function is defined as 

/>OC 

Y(s) = / dxx s - x e- x . 
Jo 

Now consider the change of variables x = X n t 

poo 

T{s) =X s n dte , - 1 e- tXn , 
Jo 

which can be inverted to give Ca(s), the zeta function of the 3-manifold 

Y-L = JL dtt^Y^ 
U(s) = ri7)J dtt Y ^)- 

In this last formula, 

Y A (t) = Y e ~ tXn > 

n 

is the trace of the heat kernel K(t, x, x'), 

d 

-— K(t,x,x') =A x K(t,x,x') 

X, X') - Y ^K^Kix') 
n 

Y A {t) = [ d s x^gK{t,x,x). 

Where the subscript x means that the operator A do not act on the variable x' , 
g and Vm ls the determinant of the metric and the volume of the 3-manifold, 
respectively. The parameter t above is not to be confused with time. 

3 Quantum Field Theory in The Static Universe 

R x H 3 

The hyperbolic space sections H 3 , can be realized as a surface 

(x - x'f + (y- y'f + (z- z'f - (w - w'f - -a 2 , (1) 

embedded in a Minkowski 4-space 

dl 2 = dx 2 + dy 2 + dz 2 - dw 2 . 

As this space is homogeneous, we explicitly write the origin of coordinates 
(x',y',z',w'). It can easily be seen that its isometry group is the proper, 
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orthocronous Lorentz group SO' (1,3), which is isomorphic to PSL(2,C) — 
SL(2, C)/{±J}. With the constraint of Eq. (QJ on the line element we obtain 



dl = dx + dy 2 + dz~ 
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{{x — x') dx + (y — y') dy + (z — z') dz) 
{x - x'f + (y - y'f + {z - z'f + a 2 ' 
ds 2 = -dt 2 + dl 2 = g{x 1 x') lLV dx^dx v , (2) 

where we interchangeably write (a: , x 1 , x 2 , x 2 ) < — > (t,x,y,z). Both connec- 
tions, V x and V x ', compatible with the metric of Eq. 0, can be defined 
through 

V M 5(x,x') Q /3 = , (3) 
V^g(x,x') af} = . (4) 

The expression in Eq. is the popular Robertson- Walker line element, which 
written in the Lobatchevsky form reads 

ds 2 = -dt 2 + a 2 [d X 2 +sinh 2 X (de 2 + sin 2 Ode/) 2 )] , (5) 

with 

sinh z x = 



{x-x'f + (y-y'f + {z-z'f 



As is well known, the EQ for the homogeneous and isotropic space sections 
in Eq. ©, with a = a(t), reduce to the Friedmann-Lemaitre equations 



d\ 2 1 8ttG A 



2 3 P 3 ' 
2[ a -) + ( a -) 2 -\ = -S,G P + A 



where the right-hand side comes from the classical energy-momentum source for 
the geometry, T^ v = (p + pju^v^ +pg^", plus the cosmological constant term 

We assume that the universe was radiation dominated near the Planck era, 
hence p = p/3, and obtain the following static solution 



/ 3 
2|A| 



_ A 

ds 2 = -dt 2 + a 2 [d X 2 +sinh 2 X {d6 2 + sin 2 Ode/} 2 )] , (6) 

where the cosmological constant is negative. 

We now wish to evaluate the vacuum expectation value of the energy density 
for the case of a universe consisting of a classical radiation fluid, a cosmological 
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constant, and a non-interacting quantum scalar field <p. The solution of EQ is 
given in Eq. J^J, where the quantum back reaction is disregarded. We use the 
point splitting method in the universal covering space R x H 3 , for which the 
propagator is exact. The point splitting method was constructed to obtain the 
renormalized (finite) expectation values for the quantum mechanical operators. 
It is based on the Schwinger formalism |.'U| . and was developed in the context 
of curved space by DeWitt Further details are contained in the articles by 
Christensen [HSj, [SZj- For a review, see |30] . 
Metric variations in the scalar action 




with conformal coupling £ = 1/6, give the classical energy- momentum tensor 

2 1 1 

+ ^g^<p\3<P + ^G^<j> 2 - -m 2 g^ 2 , (7) 

where is the Einstein tensor. As expected for massless fields, it can be 
verified that the trace of the above tensor is identically zero if m = 0. Variations 
with respect to <j> result in the curved space generalization of the Klein-Gordon 
equation, 

□0-^0-m 2 (/> = O . (8) 

The renormalized energy-momentum tensor involves field products at the 
same spacetime point. Thus the idea is to calculate the averaged products at 
separate points, x and x' , taking the limit x' — ► x in the end. 

(0\T^(x)\0) = lim T(x,x')„ v , (9) 

X 1 — >X 

with 



~ (V M V„ + V„V„) + ^g(x)^ (□ + □') 



i 



l 



+ 12 1 ~ l R ( x )9{x)»u 
GW(x,x') , 



1 



m 2 g(x)^ 



(10) 



where the covariant derivatives are defined in Eqs. <EJ and <@J, and G^ is the 
Hadamard function, which is the expectation value of the anti-commutator of 
<f){x) and (fi(x') (see below). We stress that the quantity T(x,x') flL , only makes 
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sense after the limit in Eq. is taken. The geometric quantities such as the 
metric and the curvature are regarded as classical entities. g(x) fiu = g(x,x' — 
0)^ is obtained from the line element in Eq. (J^J. 

The causal Green function, or Feynman propagator, is obtained as 

G(x,x') = i(Q\T<t>(x)4>{x')\0) , 

where T is the time-ordering operator. Taking its real and imaginary parts, 

G(x,x') = G s (x,x') + l -G^{x,x') , (11) 

we get for the Hadamard function 

G (1) (x,x') = (0\{<p{x),<j>(x')}\0) = 2lmG(x,x') . 

4 The Feynman Propagator and the Casimir En- 
ergy Density in R x j\4 

Green functions, as any other function defined in the spatially compact space- 
time iixM, must have the same periodicities of the manifold M. itself. One way 
of imposing this periodicity is by determining the spectrum of the Laplacian, 
which can only be done numerically. 

Another method imposes the periodicity by brute force, 

Im{x) = f(l x ) ■ 

The above expression is named the Poincare series, and when it converges, it 
defines functions Jm on the manifold M.. 
We define the operator 

F(x,x')=F(x)/^S(x-x') , (12) 

where F(x) = □ — R/6 — to 2 , and introduce an auxiliary evolution parameter s 
and a complete orthonormal set of states \x), such that 

G(x,x') = (x\G\x') , 
F(x,x') = (x\F\x') , 

G = i e- lsp ds . (13) 
Jo 

This last equation implies that G = (F — iO)^ 1 , hence the causal Green function 
becomes 

/•OO 

G(x,x')=i ds(x\exp(-isF)\x') , (14) 
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and the matrix element (x\ exp(— isF)\x') = K(s,x,x') satisfies a Schrodinger 
type equation, 

i^-K(s, x, x') = (u - ^ - m 2 ^ K(s, x, x') , 

where in analogy with section 12 K is named the Schrodinger kernel. 

Assuming that K depends only on the total geodesic distance — (t — t') 2 + 
a 2 x 2 , with the spatial part a 2 % 2 derived from Eq. JjjJ, the above equation can 
be solved, and we get 

, -iX e^{im 2 s + i[{t-t') 2 - a 2 x 2 ]/^s} 

A(s,x,x) = — — . (15) 

smhx (47rs) 2 

Substituting this solution for the integrand in Eq. 11411 . gives for the Feynman 
propagator 



H[ 2) (m^(t-t' f-a 2 X 2 
G(x, x 1 ) = -^-^ \ , (16) 



(t - t'Y - a 2 X 2 

(2) 

where Hi is the Hankel function of the second kind and order one. 
The Hadamard function can be obtained from Eqs. 11 (ill and lilt . 



= V (17) 



Ki ( mJ-(t-t'Y+a 2 x 2 



2ir 2 sinh x 



yj- (* - t' f + a2 x 2 



where K\ is the modified Bessel function of the second kind and order one. The 
massless limit m = can immediately be checked: 

X f 1 



G(I)(l '^ 0= 2Amh( X ){ 



2tt 2 sinh (x) I -{t - t') 2 + a 2 X 2 
Remembering that for a — > oo 



sinh x — « 1 y{% x') 2 + (y — y') 2 + (z — z') 2 — > x i 
the well known Minkowski result is recovered for the massive and massless cases, 



G (1) (x,x') = 



K x ( md - (t - t'f + r 



,■2 



^1 (i - t') 2 + r 2 
G(1) (X ' X ' )m= ° = 2^ { -(t-i') 2 + r 2 } ' 



where r is the geodesic distance in the spatial euclidean section. 
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Substituting Eq. II17H and the covariant derivatives © and J3J into Eq. 
IjlOJl . we obtain T{x,x') tiv . 

The Klein-Gordon equation remains unchanged under isometries, 



□ - | - ) , - 



□ -|-m a j i.-' 



where £^ is the Lie derivative with respect to the Killing vector £ that gener- 
ates the isometry, hence summations in the Green functions over the discrete 
elements of the group V is well defined. 

In M. = H 3 /F, a summation over the infinite number of geodesies connecting 
x and x' is obtained by the action of the elements 7 € T, which are the generators 
g% and their products (except the identity), on x' . Since Y is isomorphic to 
tti(M) each geodesic linking x and x' in M is lifted to a unique geodesic linking 
x and 72/ in if 3 . Thus from Eq. © we get 



pc = (0\T(x)^\0) M u»u» 

= u»u v lim ^ T (1,7a/) . (18) 

The infinite summation occurs because the spacetime R x M is static, so there 
has been enough time for the quantum interaction of the scalar field with the 
geometry to travel any distance. Since we know the universe is expanding, 
the infinite summation is not physically valid. The presence of the mass term, 
however, naturally introduces a cutoff. 

In Eq. II18J1 the subscript 7^1 means that the direct path is not to be 
taken into account. It can be shown that this procedure is equivalent to a 
renormalization of the comsological constant, see for instance 26j. 



5 Casimir Density p c in a Few Universes 

There are some reasons that turn smaller Universes more interesting. We de- 
scribe some spatially compact universes, with increasing volumes, in subsections 

Si 

The values of pc shown for each manifold were taken at points (8, ip) on the 
surface of a sphere inside its fundamental region. For all of them the radius of 
the sphere is the same, d = a\ = 0.390035. ..a, where d is the geodesic distance. 
Our result is displayed in Figures [3] Q] EH and |H1 for a scalar field with mass 
m = 0.4, and a metric scale factor a = 10. Angles 8 and <p correspond to the 
co-latitude and longitude, so the lines 8 = and 8 — ir correspond to the poles 
of the chosen sphere. For each manifold we also write explicitly the radius of 
the inscribed sphere Rmradius- 

The description that follows applies to all subsections I5.1I5TH The gt ma- 
trices that generate T were obtained with the computer program SnapPea 
|38j . To yield a numerical result, the infinite summation l|18ll has to be trun- 
cated. Recall that this summation occurs in the covering space H 3 . We halted 
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the summation each time the action of the generators gi and their products 
on the origin (a; 1 , a; 2 , x 3 ) = (0,0,0), reached a geodesic distance bigger than 
d — ax = 5. 29834. ..a. Care was taken so that no point was summed more than 
once. In other words, the summation l|T8l) . which yelds pc, was truncated when 
the interior of the hyperbolic sphere of geodesic radius d = a\ = 5. 29834. ..a, 
was covered with replicas of the fundamental region. We checked that additional 
contributions were unimportant for the evaluation of pc- 

5.1 Weeks Manifold 

This manifold was discovered independently by Weeks [2^ and Matveev-Fomenko 
|4H] . and is the manifold with the smallest volume (in units of a 3 ) known, 
V = 0.942707.. .a 3 . Its fundamental region is an 18-face polyhedron, shown 
in Figure^ The radius of the inscribed sphere is Rinradius = 0.519162. ..a. 

The vacuum expectation value of the 00-component of the energy- momentum 
tensor, pc — T iiV u^u v ', as seen by a comoving observer, is shown in Figure El 



Figure 1: Fundamental region for Weeks manifold. 




5.2 Thurston Manifold 

It was discovered by the field medalist William Thurston ^TJ . This manifold 
possesses a fundamental region of 16 faces, its volume is V = 0.981369. ..a 3 
(Figure El, and R lnr adtus = 0.535437.. .a. 

Figure 31 shows the value of pc — T^u^u" — Too, as seen by a comoving 
observer. 




Figure 3: Fundamental region for Thurston manifold. 




Figure 4: pc for Thurston universe. 
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5.3 Best Manifold 

This manifold was discovered as a by-product of a study of finite subgroups of 
SO(l, 3) by a geometrical aproach |32|. Its fundamental region is an icosahedron 
with V = 4.686034.. .a 3 and R lnra dzu S = 0.868298. ..a, shown in Figured 

The vacuum expectation value of pc = T^u^u" ', as seen by a comoving 
observer, is shown in Figure El 



Figure 5: Fundamental region for Best manifold. 




Figure 6: pc for Best universe. 



5.4 Seifert- Weber Manifold 

For this manifold, which was discovered by Weber and Seifert H2], V = 11. 199065. ..a 3 , 
Rinradius = 0.996384. . .a, and the fundamental region is a dodecahedron (Figure 
0). 

Figure |H1 shows the value of pc — T^u^u" ', as seen by a comoving observer. 
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Figure 7: Fundamental region for Seifert-Weber manifold. 



-O.OOUBBE-E; 



f-0,001170E-5. 



0.0011BBE-5 



0,0011 70E-5, 




Figure 8: pc for Seifert-Weber universe. 



6 Conclusions 

We explicitly evaluated the distribution of the vacuum energy density of a con- 
formally coupled massive scalar field, for static universes with compact spatial 
sections of negative curvature and increasing volume: Weeks, Thurston, Best, 
and Seifert-Weber manifolds. As a specific example, we chose m = 0.4 for 
the mass of the scalar field, and a = 10 for the radius of curvature. The val- 
ues of the Casimir energy density pc on a sphere of proper (geodesic) radius 
d = 3.90035.. inside the fundamental polyhedron for each of these manifolds are 
shown in Figures 0] EH and|H| In all these cases it can be seen that there is a 
spontaneous generation of low multipolar components. As expected, the effect 
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becomes weaker for increasing volume universes. 
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